pub fn new_birthday_probability(n: u32) -> f64 {
    1.0 - ((366 - n)..366).map(|i| i as f64 / 365.0).product::<f64>()
}
